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FOREWORD 

This report was prepared under Contract NAS8-29270 for the NASA 
Marshall Space Flight Center, and It represents a partial fulfillment 
of the requirements of that contract. The work was conducted under 
the technical cognizant of Dr. Terry Greenwood and Mr. David Seymour 
of NASA/MSFC and It was performed In cooporatlon with B. Audeh, 

R. E. Carter, and S. Smith of LMSC-HRtC. 
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Section 1 
INTRODUCTION 

The base region of the Space Shuttle Orblter receives significant thermal 
radiation from the plumes of the Space Shuttle Main Engines (SSME). This re- 
port describes the methods used In predicting this radiation and some of the 
results obtained for use as a portion of the 1978 thermal environment update. 
The complete base thermal radiation environment Is presented In Refs. 1-4. 

The procedure followed Is similar to that described In Ref. 5 for the 
1976 environment, but the plume property models have been Improved. The method 
uses plume property data at four altitudes with the NASA/MSFC gaseous plume 
radiation program (GASRAD) to predict radiation to representative surface 
locations in the base region. These predictions are used to develop (1) a 
sea level plume radiation model using a number of geometric surfaces with 
specified emissive powers and (2) an rltltude adjustment procedure to be used 
to adjust the sea level radiation predictions to the proper altitude value. 

The plume model Is used with a radiation view factor (RAVFAC) program to pre- 
dict sea level radiation at all the specified body points In the base region 
of the orblter, external tank (ET), and solid rocket boosters (SRB). These 
RAVFAC results are reviewed for consistency In comparison with expected trends 
and available GASRAD results, and they are adjusted based on the comparative 
GASRAD values when It appears that the RAVFAC plume model has given erroneous 
results. The resulting corrected sea level predictions are then published 
along with altitude adjustment factors In the thermal radiation environment 
for the SSME plumes. 

Radiation for the SRB plumes is predicted separately using RAVFAC models 
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at several altitudes developed in Ref. 6. No interaction between the SRB and 
SSME plumes is considered during development of the RAVFAC models, but opaque 
surfaces of the plume models do cause shading of other plumes in the RAVFAC 
predictions. An evaluation of the radiation potential of the SRB/SSME plume 
interaction region indicated it was not a significant radiation source. The 
evaluation and estimates of recirculating gas radiation are discussed along 
with the SSME plume predictions In Sections 2 and 3. 

The procedure of using the RAVFAC program has been used primarily to reduce 
required computer time since it is much faster than the GASRAD program. Al- 
though a significant amount of time is required to prepare the RAVFAC plume 
model and provide GASRAD predictions to check and adjust the RAVFAC results, 
an overall time saving can be realized when computing radiation to the several 
hundred points required for the base region environment of the orbiter, ET, 
and SRB. 

The plume predictions and some of the RAVFAC results presented here 
were not part of this work, but are taken from Refs. 1 thru 4 and 7 thru 9. 

The plume properties in the form used as input to the GASRAD program are pre- 
sented in this report along with a brief review of the methods used for the 
plume predictions to provide sufficient data to substantiate the results. 

The report also contains a brief description of the GASRAD program used and 
the predictions obtained, the RAVFAC plume model development, and a discussion 
of the recommended method of computing the contribution of the SSME plumes to 
the Space Shuttle base thermal environment. 
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Section 2 

SSME PLUME DESCRIPTION 

Two methods were used for plume predictions. The low altitude predictions 
(0, 20, and 40 kft) considered the atmospheric mixing and resulting combustion 
with the fuel rich plumes, but each plume was treated as an Individual axl symmetric 
plume with no interaction with the adjacent plumes. The high altitude pre- 
diction for vacuum conditions had no free stream effects considered, but did 
account for the interaction between adjacent engines with approximations for 
the complex shock structure. Complete listings of the plume properties used 
are in the Appendix. 

2.1 LOW ALTITUDE PREDICTIONS 

The low altitude plume property predictions were taken from Refs. 7-9. 

The viscous mixing analysis used a turbulent kinetic energy model, and the 
calculation was performed using the LAMP program (Ref. 10). Finite rate chemistry 
was used for the mixing analysis In all cases. The properties In the Invlscld 
portions of the plume were predicted using a method-of-characterlstlcs (M0C) 
program (Ref. 11). The mixing model was run starting with constant velocity and 
pressure for the engine flow.and a free-stream velocity of 250 ft/sec was used 
at all altitudes to represent the low velocity in the base region. The pressure 
and velocity of the plume for the mixing model were taken to be the exit wall 
boundary conditions predicted by the M0C program. The results of the mixing 
model and M0C engine plume were matched by hand. The Mach disk was located 
in the sea level plume based on photographs of test firings and In the 20 kft 
plume using the triple-point matching criteria. 
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The SSME engine arrangement shown In Fig. 2.1 shows the relatively tight 
spacing of the engines, particularly engines 2 and 3. This close spacing between 
the plumes will cause plume Interference, but this was neglected In developing 
the low altitude plumes since there were no procedures to predict the 1; ~v.r- 
ference and nonsymmetrlc mixing process. The Interference between plumes is 
expected to Impede the mixing process so that the assumption of no Interference 
would be conservative for mixing effects, and over the altitude range considered, 

there should be no significant compression due to Invlscld flow effects. 

Isotherms for the plume predictions are Illustrated In Figs. 2.2 and 2.3 

and tabulated In the Appendix. The plumes are described as "revised" to distin- 
guish them from the TKE plumes used In the 1976 prediction (Ref. 5). The sea 
level plume has a very high temperature region caused by a Mach disk at 52 Inches. 
The predicted Mach disk at 20 kft occurs at 135 Inches, but Is too small for 
the maximum temperature to be shown In Fig. 2, At 40 kft, the Mach disk Is absent, 
but the compression region on the center line has an initial peak at 160 Inches. 

The current low altitude SSME plumes have high temperature mixing regions 
which grow more slowly initially than those used in the previous environment (Ref. 5), 

but the high temperature compression region following the Mach disk In the sea 
level plume persists farther downstream. The combined effect In the sea level 
plume produces radiation which Is 1 to 20 percent lower than the previous plume 
model at representative base points shown In Table 2.1. At the two higher 
altitudes where the mixing layer is the predominate radiation source, the 
reductions for the new plume model are greater, ranging from 21 to 47 percent. 

2.2 VACUUM PLUME 

The vacuum plume was the same one used in Ref. 5 except that it was revised 
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100 inch diameter cylinders to 
illustrate relative plume positions 
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Fig. 2.1 SSME flight null engine arrangement 
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TABLE 2.1 

COMPARISON Or RADIATION PREDICTIONS FOR 
THE 1976 (TKE1 ) AND 1970 (TKE2) PLUME MODELS 



■■ 

Incident radiation - Rtu/ft ? -sec 

(ratio to TKE 1 ) 

Location 

Bjjl 

Sea Level 

20 kft 

40 kft 

TKE 1 

T KF 2 

TKE 1 

TKE 2 

TKE 1 

TKE 2 

Orbiter heat shield ctr 

905 

2.80 

2.55 

(0.84) 

3.14 

2.03 

(0.65) 

1.79 

1.26 

(0.70) 

RCS Pod aft surface 

837 

13.1 

12.1 

(0.92) 

8.60 

5.87 

(0.68) 

4.15 

3.05 

(0.73) 

Body Flap upr ctr TE 

233 

12.7 

11.6 

(0.91) 

6 . 53 

3.99 

(0.61) 

2.84 

1 .89 
(0.67) 

Fin TE 40''. span 

640 

14.2 

13.2 

(0.93) 

9.57 

6.71 

(0.69) 

5.00 

3.61 

(0.72) 

SSME 1 lip - aft 

7580 

29.8 

25.7 

(0.86) 

15.7 

9.61 

(0.61) 

6.76 

5.03 

(0.74) 

- top 

7480 

11.2 

10.7 

(0.96) 

6.20 

3.26 

(0.53) 

2.36 

1.45 

(0.61) 

LT base center 

8000 

2.88 

2.84 

(0.99) 

2.52 

1 .88 
(0.75) 

1.35 

1.06 

(0.79) 


Notes: 1. Plume model TKF 1 was used for the 1976 thermal environment 

described in Ref. 5. 

2. Plume model TKE 2 is the current model with isotherms shown 
in Figs. 2.2 and 2.3 and properties tabulated in the Appendix. 
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to provide property values at 5 Inch radial Intervals Instead of 20 Inch Inter- 
vals. The portion of the plume aft of the exit plane Is pieced together using 
four regions with three prediction techniques. The prediction techniques were 
tried for the Saturn S- II stage and checked against flight measurements to 
validate the prediction method. 

Because of the geometry restrictions required to simplify the construction 
of the vacuum plume. It was necessary to use an equivalent symmetrical engine 
cluster In which all engines are canted an equal amount from a cluster center- 
line. The cluster centerline chosen is pitched up 12.024 degrees relative to 
the orbiter X Q axis, and this centerline represents the Z axis of flow field coor- 
dinate system. Each engine Is canted 4 degrees away from this centerline and 
equally spaced around the centerline. The resulting orbiter coordinates for the 
engine exits (using a 157 Inch engine length) are compared to the actual lo- 
cations below. Using a symmetrical cluster allows the plume to be described 
by the properties In one sector as shown In Fig. 2.4. 


Engine 1 Engines 2 or 3 



x o 

Y o 

Z o 

x o 

Y o 

Z o 

Actual 

1595.9 

0 

486.3 

1622.8 

1 53 

369.9 

Equivalent Cluster 

1597.2 

0 

485.5 

1620.7 

t64.9 

375.5 

Difference 

+1.3 

0 

-0.8 

-2.1 

+11.9 

+5.6 


The procedure used In the present prediction was described In Ref. 14. 
The plume for the first 200 Inches downstream of the nozzle exit plane Is 
assembled as shown In Fig. 2.4 from an axlsymmetrlc vacuum plume about an In- 
dividual engine axis, a central core region resulting from the Impingement 
of all three plumes at the center, and a diamond shape Intention region be- 
tween adjacent plumes which Intersects the central core on the Inboard end 
and extends to an approximate plume boundary on the outboard end. The shocks 
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Fig. 2 A Composite cross-section of the SSMF flow field geometry 
for the vacuum plume 
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from the impingement regions cross the individual engine axes about 200 to 250 
inches downstream of the exit. At this point, the approximation of the 3-0 
plume was abandoned and an axisymmetrlc prediction was used from 200 to 595 
inches downstream of the nozzle exit. 

The free plume portion of the flow field was predicted using an M0C 
program (Ref. 11) run for near vacuum conditions. This plume was also used 
as input to the PL IMP program (Ref. 15) to predict pressures and temperatures 
on a surface representing the impingement plane between engines. The impinge- 
ment pressures were predicted using a modified Newtonian procedure and a 
corresponding gas temperature was predicted as that downstream of an oblique 
shock which would give the predicted pressure. If an oblique shock solution 
was not possible, the shock angle was fixed at 90 degrees relative to the flow 
and the downstream entropy and predicted pressures were used to define the 
temperature. 

The data (flow angle, Mach number, and entropy) predicted on the Impinge- 
ment plane were also used with a 2-D M0C solution to estimate the properties 
in the central core of the cluster caused by inward flow from the plumes and 
plume impingement regions. The location for the origin of the shock nea the 
cluster centerline used in the M0C solution was estimated using a modified 
stand-off distance calculation. The shock shape develops., jy the M0C solution 
was estimated using a modified stand-off distance calculation. The shock shape 
devtioped by the M0- solution using the Impingement plane flow as Input is 
shown in Fig. 2.5. The curved shock was Incompatible with tii<_ requirements 
of the DTRANS program (Ref. 16) used to assemble the flow field components, so 
a shock angle of 15.83 degrees was used as shown in Fig. 2.5. 

i he three flow field pieces (free plume, Impingement region between 
adjacent engines, and central core region) were assembled by the DTRANS 
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program into a 3-D flow field for the sector shown in Fig. 2.4 The flow 
field In the sector consists of a three-dimensional array of data arranged in 
planes normal to the flow field Z axis (cluster centerline) with radial 
values at each Z-cut at particular values of n- The radial values are 
measured from an axis originating at the nozzle exit centerline, but parallel 
to the cluster centerline. The shock locations separating the regions are 
noted In the resultant flow field data. Since the DTRANS program uses a single 
shock angle to describe both the central core and Impingement zone between 
engines, the 15.83 degree angle approximating the central core was used for 
both regions. 

Due to the geometry of the 3-D flow field description It was difficult 
to describe the central core boundary as a cone as It had original ’y been 
visualized although the true physical shape Is undefined. The two boundaries 
of the central core shown in Fig. 2.4 are a radial line and the impingement 
region boundary w ;, 1ch bo:h terminate at positions defined by the central core 
shock. The maximum depth of the impingement region on a line between the 
free plume axes is also determined by the 15.83 degree shock, but since this 
region starts nearer the engine exit than the central core, the distance to 
the shock is larger at a particular downstream location than the central core 
radius. The outer end of impingement zone between engines tapers to a point 
near the free plume boundary. 

Properties at each axial position in the central core region are constant 
radially approximating the results of the MOC solution in the same area. 
Properties at the central plane of the Impingement zone predicted by the PUMP 
program were assumed constant normal to the impingement plane. As mentioned 
previously, the 3-D flow field terminates at Z=200 Inches and the MOC prediction 
of the flow along the impingement plane and into the central core flow Is 
used alone beyond this point. 
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Typical property variations for both the 3-D flow and the MOC prediction 
are shown In Figs. 2.6 and 2.7. The solid line Indicates properties In a 
plane through an engine centerline. The constant property central core ends 
with the shock boundary and Is followed by the free plume. In a plane between 
engines, the constant property central core Is followed by the properties on 
the Impingement plane between engines represented by the dashed line. The radial 
property variation corresponding to the MOC solution used downstream of Z=?00 
inches Is represented by the broken line. It Is Interesting that the free 
plume has higher pressures at this point although the Impingement regions 
have generally higher temperatures. 

A representation of the base gas between the heat shield and engine 
exit plane was chosen to correspond approximately to current estimates of 
the recirculated flow properties on and above the heat shield. This repre- 
sentation was added to the 3-D plume flow field. The model Is based on esti- 
mated pressures using Isentroplc flow with an Isentroplc exponent of 1.3, 
a stagnation temperature of 3400R, and a stagnation pressure of 0.033 psla 
The model chosen was Intended to be more representative of the region around 
the upper engine since the base flap causes some differences In the flow In the 
region between the lower engines. The base gas model described In the following 
table uses the geomerty shown In Fig. 2.4 with n and R measured relative to the 


BASE GAS MODEL 


n 

deg. 

Max. 

radius 

inches 

Heatshield 
Z = -123 
P T 

psla R 

End RCS Pod 
Z = -38 
P T 

psla R 

Exit Plane 
Z = -0.01 
P T 

psla R 

0 

0-75.2 

0.033 

3400 

0.021 

3063 

0.011 

2639 

50 

0-69.3 

0.033 

3400 

0.021 

3063 

0.011 

2639 

78 

0-97.3 

0.02 

3029 

0.01 

2581 

0 

0 

110 

0-59.0 

0.02 

3029 

0.01 

2581 

0 

0 
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Engine centerline at a radius of 



Fig. 2.6 Distribution of static pressure in the transition 
between 3-D and 2-D plume models at Z=200 inches. 
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Engine centerline at a radius of 



Fig. 2.7 Distribution of static temperature in the transition 
between 3-D and 2-D plume models at Z=200 Inches. 


16 



F^> IS/I TF-ICH I r'-J cz. 


RTR 014-7 


upper nozzle exit centerline. The Inner region represented by n = 0 to 50 
degrees represents conditions within the cluster while the region represented 
by n = 78 to 110 degrees represents the outer base region. The radiation 
program will treat the zone between these regions as a transition zone with 
properties at any Z varying linearly with n. 

A complete listing of the vacuum plume Is In the appendix. 


2.3 OTHER PLUME EFFECTS 

In addition to the low altitude and vacuum plumes, consideration was 
given to the effects of radiation from the SSME/SRB plume interaction regions 
and the reversed gas from the SSME/SRB plume interactions. The radiation from 
these sources was not generally significant, but the effect of the reversed 
gas radiation was considered in assigning altitude correction factors for the 
SRB. 

The SRB and SSME plumes intersect in the altitude range from 60 kft 
to 140 kft as shewn by typical intersection curves in Fig. 2.8. The analysis 
of radiation from the intersection was performed at 100 kft. Although the 
surface of the intersection is not flat and has no symmetry, it was necessary 
to approximate it with a plane surface as shown in Fig. 2.9a. Properties 
on either side of the XZ plane were assumed to be symmetric about the X 
axis. The angular orientation of the plane was varied until the impingement 
pressures predicted on both sides of the plane were approximately equal. 

The ZMAX position indicated In Fig. 2.9b was used to approximate the Indi- 
cated shock position relative to the plane, and the properties were assumed 
to be constant between the impingement plane and the upper bounding line 
between ZKAX and YMAX. The estimated gas properties are listed in the 
appendix. 


1'/ 


\ 


SRB/SRB 



Fig. 2.8 Intersection boundaries for the SSME 
and SRB plumes (Ret. 23) 
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= X tan a 



a) General Arrangement 


+Z 


YMAX 

SSME Side 


1 


/ 


+ Y 




_ SSME side 
ZMAX- X tan 6.75° 

SSME Shock 


♦ 

ZMAX 


& 


SRB Side 


-Y 


c SRB Shock 

-SRB side 

ZMAX= -X tan 8.25° 


-Z 


b) Specification of the Impingement Region Shape 


Fiq. 2.9 Arrangement for describing the SSME/SRB interaction region 
properties in the Gaseous Radiation Program. 
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Cyl i rider 

Radius 

in. 

X o 

T 

R 

P 

p a fa 


Mole 

tractions 


H r 0 

CO/ 

CO 

HC1 

N? 

Sma 1 1 

54 

1613 

1560 

66 

0.03 

0.01 

0.01 

0.01 

"o?77 

(Reversed jet) 


2112 

4100 

66 

C. 33 

0.06 

0.13 

0.09 

0.06 

Large 

181 

1308 

1560 

66 

0.03 

0.01 

0.01 

0.01 

0.72 

(Separated flow) 
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Fig. 2.10 Reversed gas model for 80 kft. 
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Fig. 2.11 Reversed gas model for 120 kft. 
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Intersection of the SSME/SRB plumes also causes reversed flow in the lower 
base region. Properties and size of the reversed flow regions were conserva- 
tively estimated (Ret. IS) from analytical predictions and experimental results. 
The estimated gas volume and properties are illustrated in Figs. 2.10 and 
.11. Components of the gas property model were considered to be a small 
cylinder representing the reversed jet toward the base and a larger cylinder 

representing the separated region aft of the ET base. Temperature at the 

origin of the reversed jet (4100 R) is a prediction of the average stagna- 
tion temperature of the reversed portion of the SRB boundary layer. Tem- 
perature and pressure at the forward end is based on model measurements 
scaled to flight conditions. Gas at the origin of the reversed jet is esti- 
mated to be 66 percent SRN products and 34 percent SSME products. This gas 
is assumed to mix with sufficient air to lower the gas temperature to the 

value measured at the forward end. A linear variation of properties is 

assumed in the reversed jet and the properties in the separated region are 
taken to be equal to those at the forward end of the reversed jet. 
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Section 3 

GASEOUS RADIATION PREDICTION'S 

Predictions were made to many points in the base region using the 
gaseous radiation (GASRAD) program to serve as a basis for adjusting pre- 
dictions made by the radiation view factor (RAVFAC) program and to provide 
data on the effects of altitude on plume radiation. The essential features 
of the program used and the results obtained are presented in this section. 

3.1 program description 

The program (Ref. 12) is similar to earlier radiation programs developed 
for NASA/MSFC, but it has been designed to handle specific requirements 
of the Space Shuttle Main Engine (SSME) plumes. The program has two sets of 
subroutines. The first set prepares property data along lines-of-sight 
through the plumes while the second set performs radiation calculations 
using a band model along the lines-of-sight and integrates the result. The 
first, set of subroutines were changed to allow additional flow field input 
formats, reduce the number of radiating species considered, simplify input, 
and reduce run time. The second set of subroutines were modified to provide 
a separate subroutine for water vapor radiation with reduced output options 
and to provide for performing calculations with band model parameters 
averaged over large spectral intervals. Predictions of the SSME/SRB com- 
bined plume effects used the multiconstituent radiation prediction capability 
and special modifications of the flow field handling subroutines to handle 
the unusual data formats. 

The calculation procedure and spectral averaging method used are described 
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in Ref. 12. The radiation equation for the incident radiant flux at a point 
which the program integrates is 

AN° cos 0 sinQ A4>A0 

i.i 

6 . 4 . <>' . 0 

i ' i i 

where 0 and 4* are the hemispherical elevation and azimuth angles, w is the 
wavenumber and L is radius of a hemisphere over the point (the upper limit 
of the distance S along the lines of sight). The transmissivity x is determined 
using the statistical band model described in Ref. 17 as the single-line- 
group (SLG) model which uses the Curtis-Godson approximation to include the 
effects of inhomogeneous gases. The spectral integration is performed using 

integrated values of the Planck function (AN°) to eliminate errors due to 

u\ 

representing it by a central value when large spectral intervals are used. 

Tabulated values of absorption coefficient (k^) and line density (l/d () ) 
for water vapor used in the band model were obtained from two sources. Most 
of the data wore obtained from Ref. 17, but in the 2500 to 4500 cm -1 spectral 
range, later data compiled by Dr. Stephen Young (Ref. 13) was used for SSME 
predictions. Dr. Young's data have improved representations of the absorption 
coefficient at low temperatures while using data of Ref. 17 at high tempera- 
tures. The data are available at 25 cm -1 intervals, but they were usually 
averaged over 400 cm* 1 intervals from 1000 to 9800 cm" 1 with 

_ io 2 

k = (SVisT •C i “/ E 

- r Uj . t0 2 V 

Vd- 

W 1 
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where N 0 is the Planck function and Aw is 25 cm" 1 . The resulting values 

U) 

are given in Table 3.1 for the 1000 to 9800 cnr^range used. Evaluation of 
the results obtained with this procedure is made in the next section. 

The only radiating gas considered was water vapor and the only identi- 
fied foreign gas was nitrogen. The broadening effect of all other gases 
was lumped together, so the collision line width (cm -1 ) becomes 

Y c = p[ (0. 44(273/T) 03 + 0.09) F w + 0.09 F N 
+ 0.05(1 - F w - F N )](273/T)°- S 

where p is the pressure in atmospheres, T is the gas temperature in K, 
and F w and F^ are the water vapor and nitrogen mole fractions. 

During the spatial integration, the hemispherical elevation angle in- 
terval, A0, is held constant while the azimuth interval, A<|), can be varied to 
maintain a constant value of the parameter A6A<fsinO. This allows the spacing 
between lines in the azimuth direction to be set equal to the spacing in the 
elevation direction (selection of the values used are discussed in Section 3.2). 

The length of integration intervals along the 1 ine-of-sight (AS) is also 
varied by the program. Initially gas properties are interpolated at constant 
intervals (usually 3 to 6 inches) along each 1 ine-of-sight, but when the 
radiation computation is made, the available property values are averaged 
over larger intervals before a step in the radiation integration is made. 

This process is controlled by the input parameter TDIFF which is approximately 
equal tc the change in temperature allowed between the first and final 
temperature values in an integration interval. The interval length with the 
averaged property values is used as a step with constant properties in the 
radiation integration. 
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H20KGI 

U AVC NO. 

( 1/CM ) 
1200 . 1 
If 00. 5 
2000 . 2 
2400. 1 
2800. 1 
3200 . 3 
3600. 1 
4000. 1 
4400. 2 
4800. 5 
5200. 1 
S6CC. 2 
6000. 1 
6400. 4 
6800. 1 
7200. 1 
7600. 3 
8000. 3 
8400. 1 
8800. 1 
9200. 2 
9 6 00. 0 

H201(G ) 

WAVE NO. 

< 1/CN!) 
1200 . 2 
1600. 1 
2000 . 1 
2400. 3 
2800. 1 
3200. 6 
3600. 1 
4000. 7 
4400. 9 

4 8 CO . 2 
5200. 1 

5 6 CO . 1 
6000. 3 

64 DO . 3 
6800. 1 
7200. 1 
7600. 2 
8000. 4 
8400. 2 
8800. 1 
9200. 1 
9 6 CO • 0 


TABLE 3,1 

BAND MODEL PARAMETERS AVERAGED OVER 
40Q cm" 1 INTERVALS 

ABSORPTION COEFFICIENTS 'atm 1 cm at indicated temperature) 


TEMPER ATURES 

IDEG K ) 




300. 

60 0 . 

1000. 1550. 

?D0 0 . 

? C DP» 

7COO. 

.476-03 

2 . 232-02 

4.463-U? 5.496-02 

5.304-0? 

4.534-02 

4 . 699-02 

.463-01 

2 . 36 7-01 

1.252-01 6.9 34-02 

4 .4 06-0 7 

3.244-02 

2.458-0? 

.935-02 

3 . 06 5-02 

2.529-02 2.036-02 

1 .698-0** 

1.41 *-02 

1.237-0? 

.250-04 

1 .591-04 

4.013-04 7. P 40- 04 

1 .774-07 

2.486-03 

2.522-03 

.143-04 

2 . 160-04 

4.792-04 1. 3 3 1 - °3 

2.647-07 

*>.025-03 

4. 154-Q 7 

.510-03 

4 . 297-03 

1.126-0? 1.696-02 

2.122-C? 

2.084-02 

1 . 727-0'’ 

.241-01 

1 . 35 7-01 

•l.u?3-01 5.646-02 

3.38 6' — O'* 

2. 766- n 2 

1 . 821-0? 

.867-01 

6 . 958-02 

3.931-02 2.310-0? 

1 .574-0? 

1.228-02 

1,090-02 

.144-05 

5 . 686-05 

2.331-04 4.4t 7 -C4 

8 .348-04 

1 . 31 5-03 

1 .498-07 

.262-05 

1 . 64 4-04 

6. 071-04 1 . 777-03 

3.1 7 7 -0 7 

3.862-03 

3.458-0* 

.437-02 

1 . 377-02 

1.244-0? 8.892-03 

6.5 27-07 

5.187-03 

3. 970-07 

.127-02 

8 . 697-03 

6.071-03 4.477-03 

3.679-07 

3. 71 1-03 

2.831-0* 

.495-05 

2 . 230-05 

6.891-05 1.731-04 

4 .1 5 8-04 

0.396-04 

9.445-04 

.118-05 

1 . 462-04 

3.187-04 6.157-04 

1 .246-07 

1. 358-03 

1 . 674-03 

. 308-03 

1 . 965-03 

3.261-07 3.333-03 

3.1 69-07 

3.108-03 

2. 918-03 

.724-02 

1 . 812-02 

1.095-0? 4.956-03 

3.264-07 

2 . r 95-°3 

2. 192-07 

.293-03 

5 . 191-04 

8.334-04 7.129-04 

6.24 1-04 

5.565-04 

6. 153-04 

. 340-06 

6 .010-06 

2.121-05 7.211-05 

i .9 i n-04 

2.497-04 

3. 377-04 

.925-05 

3 . 107-05 

9.499-05 2.174-04 

3.799-04 

3.896-04 

4.434-04 

.484-05 

1 . 396-05 

6.584-05 1 .982-04 

3 .498-04 

3.179-04 

3. 054-04 

.558-06 

2.271-06 

5.721-06 2.660-05 

8.8 4 6-05 

1.006-04 

1 . 473-04 

.000 

o . no n 

0.000 0.0 CO 

1 .076-04 

1 .152-04 

1 . 196-04 

FINE 

: STRUCTURE PARAMETERS (cm' 1 ) 




TEMPER ATURES 

(DLG K ) 




300. 

600. 

1 uOL • 1500. 

2000. 

2 500. 

3000. 

.739-01 

4 . 151-01 

1.026 + 00 7. 058 + 02 

1,065+01 

3.499+01 

1 . 265+02 

.461-01 

3 . 855-01 

7.975-01 1.947+00 

4 .701+00 

1.115+01 

2.463+01 

.161-01 

2 . 99 2-01 

6.704-01 1.861+00 

4 ,8 o 3 +on 

1 . 7 8 * + Pl 

4.015+01 

.061-01 

6.08 3-01 

1.621+00 4.1 27+ 00 

9.153+00 

4.C5P+01 

1 . 939*02 

.632+00 

2 . 379+00 

3.16o*00 4.636+00 

1 ,049 + ri 

3.850+01 

1.084+0 7 

.936-01 

1 . 423 + 00 

2.263+00 2.6PR+20 

5.221+00 

1 . 752+01 

2.558+01 

• CO 9+ DO 

1 . 420+00 

1.861+00 2.650+00 

3 .36 5 + 00 

7.887+00 

1 . 369+01 

.882-31 

1 . 170+00 

1. 512 + 00 2 • ? 4 1 + n 0 

5.325+00 

1 .15 1 + 01 

1. 549+01 

.408-01 

1 . 103 + 00 

1. 756+00 4.899+00 

2.85 2 + 01 

7. r 04+01 

1 . 700+0? 

.771-01 

1 .051-01 

3. 1 71-0 1 1 . 2 30 ♦ n 0 

5.308+nn 

2.22 2+91 

8. 313+01 

.163-01 

8 . 91 9-02 

2.616-01 1.139+00 

4 .6 6 2 + 00 

1.978+01 

1 .082+0? 

.288-01 

2 . 239-01 

3.879-01 1. 557+00 

a .065+00 

3.636+01 

2. 794*02 

.897-01 

1 . 23 9-01 

8.185-01 5.961+03 

7.212+01 

2.787+02 

1 . 396 + 07 

.930-01 

8.677-02 

5.201-01 3. 387+00 

4.519+01 

1.708*02 

9. 363 + 0? 

.529-01 

6 . 90 7-02 

3.075-01 1.917+00 

1 .5 17+01 

6.210+01 

2. 771+0? 

.257-01 

3 . 22 9-02 

1.952-01 1.465+00 

7 .9 3 7 + 00 

3.556+01 

1 . 3C0 + 0? 

.339-01 

1 . 997-01 

1.U82+00 6 • 0 2 0 ♦ P G 

2.892+C1 

1.230+02 

4.115+32 

.942-01 

1 . C6 5 + 00 

2. 760+00 7.888*00 

1 .878 31 

3 . 7 04* 01 

6.068+01 

.292-01 

4 . 65 0-01 

1.274+00 3.510*00 

8.296+00 

1 .64 3+01 

2.727+01 

.161-01 

2 . 585-01 

6.794-01 l. p 31+00 

4 .57 7 + 00 

9.026+00 

1 . 433 + 01 

.603-01 

4 . 208-01 

I.U34+00 3.108+00 

8.636+00 

1 • 76 r + 01 

2.913+31 

.000 

0 . COO 

C.000 0.000 

1 .94 5 + 01 

3.P26+°1 

6.273+01 
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The geometric methods used in locating the points in the plumes at 
which properties are desired on each 1 ine-of-sight are straight forward. 

Once the point in the plume is located, linear iterpolation is performed 
parallel to the plume axis, and, in the three-dimensional plume, at constant n. 
Exceptions are made where a significant shock is present and can be located 
in the flow field. In these cases, the points used in the interpolation 
are forced to stay on the same side of the shock as the point being defined. 

In generating properties in the low altitude plumes the distance from 
a point on a 1 ine-of-sight to the axis of each of the three plumes is determined 
and the property determination is made in the closest plume. In areas where 
the plumes overlap because of the close spacing, this method essentially 
discards the portion of each plume beyond a central surface between adjacent 
plumes. 

3.2 PREDICTED RADIATION 

Radiation predictions for the SSME plumes will be described first followed 
by predictions for the SSME/SRB impingement regions and the SSME/SRB reversed 
gas. 

Predictions for the SSME plumes were made at representative points in 
the base region for the four altitude conditions considered. Predictions 
for the vacuum plume were not significantly different form the results of 
Ref. 5 which used a coarser mesh representation of the same basic flow field 
(20 inch instead of 5 inch radius increments as in the present '-acuum flow 
field). Because of the slight change in results, the vacuum plume runs 
we^e only repeated for a few points. The current results for the SSME plume 
radiation at four altitudes are presented in Table 3.2 along with the vacuum 
plume results from Ref. 5. The maximum predicted radiation in the GASRAD 
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TABLE 3.2 

RADIATION PREDICTIONS AT SELECTED LOCATIONS 


Location 

O 1 

Point 

No. 

- — — M 

Incident radiation - Btu/ft z -sec and 
(ratio to sea level) 

Sea 

Level 

20 

Kft 

. . -,j 

40 

K*t 

Vacuum 

Vacuum 
(Ref. 5) 

Orbiter base heat shield 




” -H 




Bottom center 


900 

3.53 

2.25 

1.23 

0.47 

0.48 





(0.64) 

(0.35) 

(0.13) 

(0.14) 

Center 


905 

2.35 

2.03 

1.26 

0.74 

0.73 





(0.86) 

(0.54) 

(0.31) 

(0.31) 

Bottom outboard 


909 

2.10 

1.78 

1.15 

0.04 

0.05 





(0.85) 

(0.55) 

(0.02) 

(0.02) 

Bottom No. 1 bulge 

270 

941 

2.47 

1.83 

1.01 






(0.74) 

(0.41) 



Top outboard 


957 

0.99 

0.78 

0.55 







(0.79) 

(0.56) 



Upper middle outboard 


961 

5.02 

3.23 

1.79 

0.42 

0.46 




(0.64) 

(0.36) 

(0.08) 

(0.09) 

Bottom eng mtd hs I 

270 

984 

3.45 

2.09 

1.03 

0.67 

0.67 




(0.61) 

(0.30) 

(0.19) 

(0.19) 

RCS Pod 








Lower surface 


807 

0.86 

0.38 

0.20 


0.04 





(0.44) 

(0.23) 


(0.05) 

Aft surface 


837 

12.1 

5.87 

3.05 

0.69 

0.70 





(0.48) 

(0.25) 

(0.06) 

(0.06) 

Outboard lower surface 


851 

0.45 

0.17 

0.10 


0.003 





(0.38) 

(0.22) 


(0.01) 

Inboard surface 


862 

4.27 

1.62 

0.77 

0.30 

0.32 





(0.38) 

(0.18) 

(0.07) 

(0.07) 

Top 


874 

0.030 

0.038 

0.034 







(1.27) 

(1.13) 



Top 


877 

0.105 

0.110 

0.081 


0.003 





(1.05) 

(0.77) 


(0.03) 

OMS Pod 








Base 


784 

4.71 

2.80 

1.57 


0.10 





(0.59) 

(0.33) 


(0.02) 

Nozzle shroud base 


7780 

7.41 

3.51 

1.69 


0.52 





(0.47) 

(0.23) 


(0.07) 

Body flap trailing edge 








Upper center 


233 

11.6 

3.99 

1.89 

0.93 

0.96 




(0.34) 

(0.16) 

(0.08) 

(0.08) 

Aft center 


241 

11.5 

5.28 

2.70 


0.58 




(0.46) 

(0.23) 


(0.05) 

Fin trail ing edge 








40'. span aft 


640 

13.2 

6.71 

3.61 

0.53 

0.59 




(0.51) 

(0.27) 

(0.04) 

(0.04) 

40" span side 


6499 

1.63 

0.93 

0.54 


0.072 





(0.57) 

(0.37) 


(0.04) 
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RADIATION PREDICTIONS AT SELECTED LOCATIONS 


Location 

0 1 

Point 

No. 

Wing trailing edge 
60/.' span aft 


360 

60" span top 


4699 

OMS Nozzle 
Aft facing lip 

180 

7150 

Inbd facing exterior surface 

180 

7152 

Inbd facing inner 1 ip 

0 

7201 

Outbd facing inner lip 

180 

7251 

RCS Aft Nozzles 8 and 9 
Inbd facing inner lip 

0 

9502 

Outbd facing inner lip 

0 

9652 

Inbd facing outer lip 

180 

9651 

SSME Nozzle 1 
Hat-band 2 top - outbd 

0 

7402 

Exit manifold top 

270 

7480 

Hat -band 2 top - inbd 

270 

7482 

Hat-band 4 top - inbd 

270 

7484 

Hat-band 2 aft - outbd 

0 

7502 

Exit manifold aft - outbd 

90 

7530 

Exit manifold aft - inbd 

270 

7580 

Hat-band 2 aft - inbd 

270 

7582 

Hat-band 4 aft - inbd 

270 

7584 

SSME Nozzle 2 
Hat-band 1 top 

135 

7741 

Hat-band 2 top 

135 

7742 

Hat-band 2 top 

18C 

7752 


Incident radiation - Btu/ft 2 -sec and 
(ratio to sea level) 


Kft Vacuum 




2 

.35 

0 

.67 

1 

.29 

0 

.59 

5 

.67 

0 

.48 

1 

.96 

0 

.33 

4 

.42 

0 

.42 

1.007 

1 

.40 


U. I I 
( 1 . 10 ) 
1.48 
(0.30) 

0.14 

(0.42) 

3.26 

(0.30) 

2.00 

(0.30) 

1.17 

(0.32) 

0.28 

(0.25) 

6.80 

(0.40) 

9.61 

(0.37) 

3.15 

(0.34) 

1.55 

(0.30) 

1.30 

(0.31) 

1.30 

(0.29) 

0.54 

(0.39) 


1.34 

0.38) 

0.70 

0.32) 

2.89 

0.24) 

0.91 

0.15) 

2.13 

0 . 20 ) 

).010 

;2.oo) 

1.72 

(0.19) 

0.08 

(0.80) 

0.64 

(0.13) 

0.08 

(0.24) 

1.45 

(0.14) 

0.88 

(0.13) 

0.50 

(0.14) 

0.17 

(0.15) 

4.06 

(0.24) 

5.03 

( 0 . 20 ) 

1.28 

(0.14) 

0.61 

( 0 . 12 ) 

0.89 

(0.15) 

0.64 

(0.14) 

0.27 

(0.19) 




063 0.093 

.02) (0.03) 

0.075 
(0.03) 


0.61 

(0.05) 


0.08 

(0.07) 


1.17 

(0.13) 

0.46 

(0.09) 

0.76 

(0.13) 


0.63 

(0.05) 

0.42 

(0.07) 

0.63 

(0.06) 

0.007 

(1.40) 

0.59 

(0.07) 

0.008 

(0.08) 

0.41 

(0.08) 

0.05 

(0.15) 

1.55 

(0.14) 

0.74 

( 0 . 11 ) 

0.51 

(0.14) 

0.08 

(0.07) 

0.79 

(0.05) 

2.86 

( 0 . 11 ) 

1.20 

(0.13) 

0.61 

( 0 . 12 ) 
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TABLE 3.2 (Continued) 

RADIATION PREDICTIONS AT SELECTED LOCATIONS 


Location 

O 1 

Point 

No. 

Incident radiation - Btu/ft 2 -sec and 
(ratio to sea level) 

Sea 

Level 

20 

Kft 

40 

Kft 

Vacuum 

Vacuum 
(Ref. 5) 

SRB Case/Skirt ring 








Aft facing web - top 

0 

6503 

3.70 

2.68 

1.53 


0.032 





(0.72) 

(0.41) 


(0.01) 

Inbd facing rim - side 

90 

6535 

2.27 

1.19 

0.61 


0.092 





(0.52) 

(0.27) 


(0.04) 

SRB stiffener at Xn=1733.56 








Aft facing web -top 

0 

6303 

2.86 

2.12 

1.23 


0.024 





(0.74) 

(0.43) 


(0.01) 

- side 

90 

6333 

3.11 

2.18 

1.24 






i 

(0.70) 

(0.40) 



Outside rim - top 

0 

6305 

3.33 

1.78 

0.95 


0.10 





(0.53) 

(0.29) 


(0.03) 

- side 

90 

6335 

1.49 

0.80 

0.41 







(0.54) 

: — — m 

(0.28) 




NOTES: 1) 0 is measured clockwise from the axis of the surface in the neqative 

Y 0 direction for the Orbiter. For the SRB and ET it is measured from 
the top (positive Z 0 ). 

2) Nozzle hat-bands numbered forward starting with the one next to the 
1 ip manifold. 
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results was at Point 7580 on the SSME engine exit lip at sea level (25.7 
Btu/ft 2 -sec). Many other relatively unshaded aft facing surfaces near the 
SSME cluster have radiation rates in the range of 10 to 13 Btu/ft 2 -sec, 
but shading by the SSME nozzles limits the peak radiation for the base heat 
shield to 5 Btu/ft z -sec. Peak radiance for a single line-of-sight through 
the plumes reached 47.8 Btu/ft 2 -sec-sr for the aft facing nozzle lip on 
SSME nozzle 1 (Point 7530). 

Although the effect of engine gimbaling was not covered as part of the 

thermal environment, sensitivity of the radiation rates to gimbaling was 
sampled by predictions of radiation to Point 7852 (the aft surface of hat 
band 2 on the inboard side of engine 2) with the lower engines gimballed 
2, 4, and 6 degrees toward the left from flight null. The results indicate 
an increase in radiation of about 12 percent ; er degree of yaw gimbal, so 
any flight conditions which are expected to have a significant gimbal de- 
flection for a long period of time should he evaluated to determine if the 
environment based on the flight null condition (shown in Fig. 2.1) is adequate. 

An evaluation was made of the effect of several parameters on the pre- 
dicted radiation. It included effects of the modified band model parameter 
data, the spectral integration interval, and the spatial integration inter- 
val along line-of-sight. The revised band model data in the 2500 to 4500 cm" 1 
range gives slightly lower radiation predictions as indicated in Table 3.3, 
so when it is combined with the previously used data outside this spectral 
range the result is to reduce the predicted radiation for the typical SSME 
plume conditions considered. The 400 cm“ l integration interval used for the 
prediction in Table 3.2 produces a slightly higher radiation than the 25 cm' 1 
interval in all cases. This offsets the apparent reduction in predicted 
radiation when TDIFF* a 200 R Is used. The TDIFF parameter extends 

* See page 25 for description of TDIFF. 
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fABLE 3.3 


RESULTS Of RAND M00FL PARAMETER 
AND PATH LENGTH INTEGRATION INTERVAL VARIATIONS 


Absorption 

Coefficient 

Date 

TDIFF 

Wave 
Number 
1 /cm 

ET Base 
Center 
Point 8000 

Nozzle Lip 
Norma 1 to 
Enq Ctrl 
Point 7480 

Fin TE 

Side 

Point 6499 

Ini t 

Final 

A 

SL 

^0 

kft 

40 

kft 

i 

I 

j SL 

20 

kft 

40 

kft 

SL 

20 

kft 

40 

kft 

NASA 




IBM 

2.91 

1.93 



3.28 

1.43 

1.67 

0.95 


(Note 1 ) 


















||M 


25 

2.87 

1.89 

1.06 










2500 

4500 

25 


0.72 

0.43 







Young 





25 



0.41 







(Note 2) 





1 










Combined 

1 


Wililll 



2.84 

1.88 

1.06 


3.26 

1.45 

1.63 

0.93 


(Note 3) 
















1 


1000 

9800 

25 

2.80 


1.04 


3.16 

1 .41 



0.53 


0 




2.87 

1.91 


10.8 

3.39, 

1.54 

1.65 

0.95 

0.56 


NOTES: 1) Data original provided by GD/'C (Refs. 20-22) and included in NASA SP 3080 

(Ref. 17) 

2) Data from Dr. Young (Ref. 13) using half-widths from Ref. 17 

3) Combination of Dr. Youngs data (2S00 - 4500 cm' 1 ) and NASA SP 3080 
data (1000 - 2475 cm" 1 and 4525 - 0800 cm' 1 ) 

4) All thermal t?nvironment predictions use the combined absorption coefficient 

data with 1QIFF = 200, = 1000, u f = 9800, and a«i= 400 
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the integration i r. t vval by averaging the property data which are available 
at G inch intervals. When TDIFF = 0, the averaging procedure i' - not used, so 
a step in the radiation integration is perfor ed at each 6 inch interval along 
a 1 ine-of-sight. 

All of the errors apparent in this comparative evaluation are insignifi- 
cant in comparison to other uncertainties in the prediction and use of the 
radiation environment, so t.he use of a 400 cm 1 spectral Interval and TDIrF = 20CR 
are considered satisfactory. 

Angular intervals used in integrating over the spherical angles are more 
difficult to evaluate since errors incurred are likely to vary significantly 
with point location. Therefore, the intervals were chosen to be approximately 
equal to the 6 inch integration interval along the 1 i ries-of-s i ght at a distance 
of about 200 inches. The cora length represented by the azimuth angle increment 
(,',>) is a function of the elevation angle (O) and increment (AO). At a radius 
or 1 ine-of-sight length (S), the cord (C) is 

C = S A,' sin 0 

for small increments expressed in radians. 

In making the SSME plume radiation prediction the procedure followed is 
to use a constant value of AO and vary AT with O to keep C/S constant*. The 
values chosen were AO = 2° and r /S = AO , so when O is 90‘ , AT = AO. The 
resulting cord length is 

C = 0.035S, 

so for AS = 6 inches, C is equal to AS at S = 171 inches. 

Radiation predictions for the SSME/SRB plume intersection region were made 
at only 4 points in the base region because the indicated heating rates wet ° 


* In the radiation program, C/S iscontroled by the input parameter TANGLE. 
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so low. The results are shown in Table 3.4 in comparison with the current 
predicted rates due to the SSME and SRB plumes at the same altitude. The results 
for the intersection regions are conservative (high) estimates since they neglect 
absorption by the SSME and SRB plumes. In all cases, radiation from only one 
(the closest) intersection region was predicted. This was doubled to give the 
results for two intersection regions for the two points on the vertical plane 
of symmetry. For the two other points, it was estimated that the second inter- 
section would increase the radiation by no more than 50 percent, so the single 
region prediction was increased by 50 percent for the comparison in Table 3.4. 

The integration intervals used for the SSME/SRB plume intersection regions 
were the same as used on the SSME plumes with the exception that the increment 
along the 1 ine-of- sight was reduced from 6 inches to 3 inches because of the 
small thickness of the initial high temperature region in the intersection flow- 
field (described in Section 2.3). 

Predictions were made at the ET base center and lower body flap trailing 
edge to evaluate the impact of radiation from the reversed gas from the SSME/SRB 
plume intersection. The results shown in Table 3.5 compared with the current 
environment indicate significant radiation could result from the reversed flow 
at 120 kft. Since the reversed flow radiation increases with altitude while 
most other radiation decreases with altitude, the altitude adjustment schemes* 
used do not furnish the necessary trends to define the effect. 

The lower body flap surface receives no SSME plume radiation, so the SSME 
environment assigned is intended to recognize a small contribution cf the reversed 
gas radiation. A full description of the radiation predicted would have required 
either artifically raising the lov' altitude heating rate (presently 0 up to 40 kft) 


* Recommend altitude adjustment procedures described in Section 5. 
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TABLE 3.4 

COMPARISON OF RADIATION ESTIMATES FOR 
THE SSME/SRB PLUME INTERSECTION REGIONS 
WITH THE DESIGN ENVIRONMENT AT AN ALTITUDE OF 100,000 FEET 


Location ; 

i 

0 

(Note 1 ) 

Point 

No. 

Predicted Radiation - Btu/ft*-sec 

Environment 
(Refs. 1-3) 

SSME/SRB 
Intersection 
Estimate (Note 2) 

SRB 

|| 

Fin trailing edge 

- 

640 

1.84 

2.64 

0.10 

ET base center 

- 

8000 

3.15 

0.65 

0.05 

SRB skirt aft rim 

0 

2100 

7.71 

1.02 

0.12 

SRB skirt inbd rim 

45 

2119 

L > 
GJ 

O 

0.85 

0.02 


NOTES: 1) 0 defined in Table 3.2 

2) Predictions made for the closest intersection region were doubled 
for points on the vertical plane of symmetry (pts. 640 and 8000) 
and increased 50 percent for the other two points to give an 
estimate for both regions. 


TABLE 3.5 

COMPARISON OF RADIATION ESTIMATES FOR THE REVERSED 
GAS WITH THE DESIGN ENVIRONMENT AT 80,000 AND 120,000 FEET 


Location 

Point 

No. 

Alt 

kft 

Predicted Radiation Btu/ft'-sec 

Envi ronment 
(Refs. 1-3) 

Reversed Gas 
E s t i ma te 

SRB 

SSME 

Body flap lower TE 

213 

80 

2.24 

0.04 

0.11 



120 

1.79 

0.07 

0.60 

FT base center 

8000 

80 

3.51 

0.81 

0.12 



120 

2.89 

0.45 

0.53 
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or ass i (jn i rui special altitude adjustment factors for this surface. Neither 
approach was doomed justified since the thermal protection system intended for 
the body flap lower surface is likely to he designed by re-entry heating loads. 

A similar anumient was made for the lower surface of the lower SSME nozzles 
where the thermal protection system is designed by high radiation and impinqe- 
ment heating at other locations on the interior of the enqine cluster. 

Radiation from the reversed qas to the FT base and the inner surface of the 
SKB were considered in assigning the SSME altitude adjustment factors and slight 
adjustments were made to the 160 kft heating rates in some cases to cover 
possible reversed gas radiation prior to the SRR separation. Surfaces other than 
those discussed here are not expected to receive significant radiation from 
the reversed gas. 

Integration intervals used in the reversed gas prediction were the same 
as those for the SSMI plumes except that the angular integration intervals 
were increased because of the relative uniformity of the oas and the close 
proximity of the gas to the prediction points. 
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Section 4 

VIEW FACTOR FLUME MODEL 

In order to allow rapid prediction of radiation distributions, a model 
representing the sea level plume was developed for use with a modified version 
of the RAVFAC program (Ref. 19) which sums the products of computed view fac- 
tors and emmissive powers for plume components to arrive at predictions of in- 
cident radiation. This section will describe the development of the plume model 
and present a comparison of the RAVFAC and GASRAD predictions for surfaces in 
the Space Shuttle base region. 

The modified RAVFAC program can handle either opaque or transparent shapes, 
so portions of the plume believed to have relatively high transmissivity can be 
modeled as transparent while those with low transmissivity can be specified as 
opaque. Because of the success achieved with the model of the previous plume 
predictions (Ref. 5), a similar model was chosen initially for development of 
the current plume model. However, it was difficult making the nozzle lip pre- 
dictions high enough without forcing the upper body flap point too high. A 
new shape was tried with a disk inside a cylinder to shade it from view for 
directions approaching normal to the plume axis. The resulting model is de- 
scribed in Fig. 4.1. In general, the size of the shapes decreased while 
emissive powers increased in comparison with the model in Ref. 5. In the present 
model, the disks were made transparent to save shading calculation time since 
they would primarily shade the interior of the cylinders which do not radiate. 

Ratios of RAVFAC to GASRAD predictions for the points used in developing 
the plume are shown in Table 4.1 for both the design plume and variations made 
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Distance from the exit-inches 


Shape 

Emissive 

power 

Btu/ft 2 -sec 

Distance 
from exit 
inches 

Radius 

inches 

6-elements 

Y -elements 

Disk 1 - Transparent 

150 

50 

25 

2 

8 

2 - Transparent 

120 

350 

20-50 

3 

8 

Cylinder 1 - Opaque 

60 

50-350 

20 

6 

16 

2 - Opaque 

60 

300-600 

50 

6 

16 

3 - Opaque 

40 

600-1200 

50 

6 

16 

Cone - Transparent 

12 

120-1200 

40-100 

11 

8 

Fig. 

4.1 RAVFAC model 

for the TKE 2 

sea level 

plume 
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to test sensitivity to element size. The goal was to achieve ratios of 1.0 to 
1.2 for the points used. The largest error (22.7%) was on the body flap upper 
trailing edge which seemed to be difficult to model. It was surmised that in 
the GASRAD prediction cooler portions of the lower plumes absorbed more of the 
high radiation levels from the upper plumes than could be duplicated within the 
RAVFAC model limitations. 

In designing a RAVFAC model, an attempt is made to minimize the number of 
surface elements, so that more computer memory will be available for surfaces 
receiving radiation and computation time will be reduced. At the same time, the 
elements must be properly distributed to allow a reasonably accurate repres- 
entation of shading effects. This also usually satisfies the requirement of 
having elements small enough for the necessary numerical computation requirements. 
The number of elements for the surfaces in the design plume were assigned 
considering these factors, then checked by varying the element size to deter- 
mine the sensitivity and detect errors due to improper sizing. The results of 
varying element size presented in Table 4.1 indicate that increasing the size 
(reducing the number of elements) does not produce a significant deteriation 
in the plume model accuracy. However, no reductions in the number of elements 
were considered since the number of elements in the plume model were not con- 
sidered excessive and there is always the possibility of larger errors at other 
locations. Increasing the number of elements did not result in any significant 
increases in the indicated RAVFAC plume error although the maximum error increased 
from 22.7 to 24.7 percent. The largest percentage change due to increasing the 
elements was at Point 8670 on the ET base where the RAVAC prediction decreased 
5 percent when elements in the axial direction on the cylinders (6 - elements) 
were doubled. This illustrates the relative sensitivity that heavily shaded 
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T"BLF. 4.1 

RESULTS OF THE RAVFAC PLUME STUDY 


Location 

Poi nt 
No. 

Ratio 

of RAVFAC/GASRAD predictions jj 

Design 

- 

3 - elements 
Note 1 

y-elements 
Note 2 




Double 

Halve 

Di sks 
16 

Cyl’s 

8 

Orbiter heat shield 
Bottom outboard 

909 

1 .098 

1.085 

1.119 

1.066 

1.102 

Upper middle outboard 

961 

1 .095 

1 .055 

1.128 

1.096 

1.115 

RCS Pod 







Inboard facing surface 

862 

1 .14* 

1.142 

1.124 

1 .162 

1.142 

Body flap trailing edge 
Upper center 

233 

1.227 

1.247 

1.211 

1.226 

1.227 

Fin trail inq edoe 







40" Span - aft 

640 

1 .030 

1.026 

1.046 

1.014 

1.040 

- side 

6499 

1.062 

1.061 

1.071 

1.046 

1.096 

Wing trailing edge 







60'P Span - aft 

360 

0.993 

0.999 

1.000 

1.001 

1.002 

SSME 1 







Exit manifold - top 

7480 

1.054 

1.069 

1.000 

1.055 

1.058 

- aft 9=90 

7530 

0.991 

1.009 

0.953 

0.995 

1.012 

- aft 0=270 

7580 

1.061 

1.070 

1.032 

1.059 

1.072 

FT Base 







Center 

8000 

1.053 

1.054 

1.070 

1 .058 

1.068 

108.4 inch radius - top 

8670 

1 .165 

1.100 

1.139 

1.153 

1.159 

ET/SRB Aft support 







I Upper link facina aft 

775 

1.139 

1.129 

1.129 

1.123 

1 .122 

1 SRB Skirt 







Outer lip - top 

2109 

1.110 

1.105 

1.136 

1.106 

1.107 

Base - top 

2500 

1.123 

1.111 

1.160 

1.121 

1.094 

Mean 


1 .090 

1 .084 

... 

1.088 

1.085 

1.094 


Notes: Changes described are variations to the "design" RAVFAC plume in Fip4.1. 

1 ) 3 -elements are the axial elements on cylinders and cones and the radial el- 
ements on disks. Since some surfaces have odd numbers of B-elements, the 
halving procedure followed was 

(B-elt)NFW=((B-elt)OLD+l)/?. 
with the result truncated to an integer . 

2) -y-elements are the circumferential elements. The yelements on the disks 
were increased to 16 to give 16 elements on all surfaces but the transparent 
cone. In reducina elements, the cylinders were reduced to 8 elements so all 
surface-; had 8 y-elements. 
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points can have to the element distribution, but it was not considered sufficient 
cause to increase the elements in the plume model. 

A comparison of the available RAVFAC and GASRAD predictions using the plume 
model are shown in Table 4.2. The indicated errors in the RAVFAC predictions 
are less than 20 percent for all significant heat rates with the exception of 
9 points on the SSME nozzles. 

Analysis of the errors indicated that too much emphasis had been placed on 
the disk located 50 inches from the nozzle exit which represents the beginning 
of the high temperature region caused by the plume Mach disk. This had been 
done to reduce the number of elements in the plume and the indicated error at 
the upper body flap point (Point 241). A better plume probably could have been 
achieved using two disks with lower emissive powers, but the errors caused by 
the use of one disk were not apparent in the sample problems used in developing 
the plume. 

The possibility of increasing the number of elements on the first disk was 
investigated to determine if this would improve the agreement for the nozzle 
points. Increasing the elements from 16 (2 x 8) to 48 (3 x 16) caused significant 
changes in some of the predictions with poor agreement, but the trend was not 
always toward reducing the error. 

On SSME nozzle 1 the agreement is good at the exit plane (Points 7XX0), 
but apparently to much radiation in the RAVFAC prediction comes from the near 
portion of plumes 2 and 3 while the mixing region of plume 1 contributes a 
negligible amount. Moving forward on the nozzle to hat band 2 (Points 7XX2) 
causes the GASRAD results to drop apparently due to shading of plume 1 by the 
nozzle exit manifold and aft most hat band, while the same location on the RAVFAC 
model does not have appreciable smaller view factors for the near plume of SSME 
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TABLE 4.2 

RADIATION PREDICTION COMPARISONS AT SELECTED LOCATIONS 


Location 

0 1 

Point 

No. 

Incident Radiation 
Btu/ft 2 sec 

Ratio 

RAVFAC/ 

RAVFAC 

.. _ 

GASRAD 

GASRAD 

Orbiter base heat shield 






Bottom center 


900 

3.74 

3.53 

1.06 

Center 


905 

2.63 

2.35 

1.12 

Bottom outboard 


909 

2.31 

2.10 

1.10 

Bottom No. 1 bulge 

270 

941 

2.61 

2.47 

1.06 

Top outboard 


957 

1.07 

0.99 

1.08 

Upper middle outboard 


961 

5.50 

5.02 

1.10 

Bottom eng mtd hs 1 

270 

984 

3.18 

3.45 

0.92 

1 RCS Pod 






1 Lower surface 


807 

0.77 

0.86 

0.90 

I Aft surface 


837 

13.4 

12.1 

1.11 | 

Outboard lower surface 

1 1 


851 

0.38 

0.45 

0.84 

Inboard surface 


862 

4.88 

4.27 

1.14 

Top 


874 

0.02 

0.03 

0.67 

Top 


877 

| 0.07 

0.105 

0.67 

: OMS Pod 



1 

| 

i 

'! Base 

- j 


784 

4.90 

! 4.71 

| 1.04 

Nozzle shroud base 

\ 


7780 

7.08 

7.41 

0.96 

1 

Body flap trailing edge 



1 



\ Upper center 


233 

14.2 

11.6 

1.23 

Aft center 


241 

| 

12.9 

11.5 

1.12 

:j Fin trailing edge 






i | *40* span aft 


640 

13.6 

13.2 

1.03 

j *40« span side 


6499 

1.73 

1.63 

1.06 

;! Wing trailing edge 






*60i span aft 


360 

3.49 

3.51 

0.99 

60 % span top 


4699 

1.79 

2.17 

0.82 

1 OMS Nozzle 






Aft facing lip 

180 

7150 

13.1 

11.8 

1.11 
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TABLE 4.2 (Continued) 

RADIATION PREDICTION COMPARISONS AT SELECTED LOCATIONS 


Location 

O 1 

"** " 

Point 

No. 

Incident Radiation 
Btu/ft 2 sec 

Ratio 

RAVFAC/ 

RAVFAC 

GASRAD 

GASRAD 

Inbd facing exterior surface 

180 

7152 

6.66 

5.93 

1.12 

Inbd facing inner 1 ip 

0 

7201 

11.6 

10.6 

1.09 

j Outbd facing inner lip 

180 

7251 

0 

0.005 

— 

RCS Aft Nozzles 8 and 9 






Inbd facing inner lip 

0 

9502 

9.75 

8.90 

1.10 

Outbd facing inner lip 

0 

9652 

0.01 

0.10 

0.10 

Inbd facing outer lip 

180 

9651 

5.82 

4.92 

1.18 

SSME Nozzle 1 



| 



!i Hat-band 2 top - outbd 

0 

7402 

0.31 

0.33 

0.94 

l 

*Exit manifold top 

270 

7480 

11.3 

1 

10.7 

1.05 

Hat-band 2 top - inbd 

270 

7482 

9.14 

6.62 

1.38 

M Hat-band 4 top - inbd 

270 

7484 

3.63 

3.68 

0.99 

1 j 

Hat-band 2 aft - outbd 

0 

7502 

1.59 

1.12 

1.42 

[ *Exit manifold aft - outbd 

90 

7530 

i 16.8 

i 

17.0 

0.99 

*Exit manifold aft - inbd 

270 

7580 

| 27.3 

25.7 

1.06 

; Hat-band 2 aft - inbd 

270 

7582 

11.0 

t 

9.20 

1.20 

i Hat-band 4 aft - inbd 

270 

7584 

! 5.20 

! 

5.23 

0.99 

SSME Nozzle 2 



I 

i 

i 

i 


Hat-band 1 top 

135 

7741 

7.56 

i 5.87 

1.29 

1 Hat-band 2 top 

135 

7742 

6.21 

4.42 

1.40 

Hat-band 2 top 

180 

7752 

! 1.02 

i 1.39 

1 

0.73 

:j Hat-band 2 top 

225 

7762 

1.64 

| 1.20 

1.37 

ij Hat-band 1 aft 

135 

7841 

12.6 

10.3 

1.22 

,! Hat-band 2 aft 

135 

7842 

9.28 

7.06 

1.31 

Hat-band 2 aft 

180 

7852 

3.93 

3.76 

1.05 

I* Hat-band 2 aft 

225 

7862 

4.00 

2.63 

1.52 

ET Base 






! ^Center 


8000 

2.96 

2.84 

1.04 

57.1 inch radius - top 

I 

0 

8330 

2.81 

2.77 

1.01 
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TABLE 4.2 (Continued) 

RADIATION PREDICTION COMPARISONS AT SELECTED LOCATIONS 


; 

i 

Location 

i 

O 1 

Point 

No. 

Incident Radiation 
Btu/ft 2 sec 

Ratio 

RAVFAC/ 

RAVFAC 

j 

GASRAD 

GASRAD 

\ 

57.1 inch radius - bottom 

180 

8339 

2.20 

2.07 

1.06 

*108.4 inch radius - top 

0 

8670 

1.92 

1.65 

1.16 

- bottom 

180 

8679 

0.93 

0.88 

1.06 

154.2 inch radius - top 

0 

8950 

0.13 

■ 

0.21 

0.62 

- bottom 

ET/SRB aft support 

180 

8959 

0.02 

f 

0.021 

0.95 

♦Upper link - aft 


775 

1.71 

1.50 

1.14 

Lower link - aft 
ij SRB aft support post 


825 

1.72 

1.71 

1.01 

' Outer lip 

'j SRB skirt 

30 

671 

9.23 

1 

! 

| 

6.98 

1.03 

ij *0uter lip - top 

0 

2109 

J 5.31 

4.78 

1.11 

-si de 

90 

2139 

1.40 

1.5Q 

1 i 

0.88 

♦Base web - top 

0 

2500 

3.68 

; 3.28 

1.12 

SRB base curtain 
SRB case/skirt ring 

45 

2317 

9.04 

l 

: 9.03 

1.00 

■j Aft facing web - top 

0 

6503 

j 4.14 

! 3.70 

j 1.12 

Inbd facing rim - side 
SRB stiffener at X^= 1 733 . 56 

90 

6535 

| 2.23 

2.27 

1 

i 

0.98 

Aft facing web - top 

0 

6303 

| 3.13 

j 2.86 

1.09 

- side 

90 

6333 

3.16 

1 3.11 

1 .02 

Outside rim - top 

0 

6305 

3.54 

! 3.33 

1.06 

- side 

90 

6335 

1.53 

j 1.49 

t — 

1 .03 


* Points used in developing the RAVFAC plume model. 

NOTES: 1) 0 is measured clockwise from the axis of the surface in the negative 
Y 0 direction from the Orbiter. For the SRB und ET It Is measured 
from the top (positive Z Q ). 

2) Nozzle hat-bands numbered forward starting with the one next to the 
1 ip manifold. 
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2 and 3. This results in an over-prediction by the RAVFAC model. The success 
of the RAVFAC model at the fourth hat-band (Points 7XX4) is fortuitous. 

On SSME nozzle 2, the over-prediction at hat-band 2 (Points 7XX2) is similar 
to that of nozzle 1 except for the points at 0 = 180° (Points 7X5X). At this 
point, nozzle 3 is so close to nozzle 2 that the initial disk surface of plume 3 
is shaded by nozzle 3. (Refer to Fig. 2.1 for engine position.) 
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Section 5 

RECOMMENDED THERMAL ENVIRONMENT 

The recommendations resulting from this study have been Incorporated In 
the radiation environment presented In Refs. 1 thru 3 to provide a sea level 
prediction with multiplying factors to account for altitude variations. 

The adjustments to the sea level predictions made using the RAVFAC program 
and plume model were generally minor except on the SSME nozzles as Indicated by 
the comparison with GASRAD predictions In Table 4.2. The adjustments were 
judiciously applied using comparisons at known points as a guide to points with 
similar orientations or locations. Because of the peer agrt ...ent between the 
RAVFAC and GASRAD predictions on the SSMl nozzles (Table 4.2) It was necessary 
to Increase the number of GASRAD predictions 1r. this region to assure the RAVFAC 
results could be properly adjusted. It Is recommended that the RAVFAC plume 
model be revised before It is used for further radiation studies on the SSME 
nozzles. 

Altitude adjustr.jnts were divided Into three ranges: 0 to 40 kft., 40 to 
160 kft., and over 160 kft. Double letter codes were provided for each point 
with the first letter designating a curve depicting the variation In adjustment 
factor from 0 to 40 kft and the second letter designating a heating rate c* 

160 kft and above. Between 40 and 160 kft a linear variation with altitude Is 
assumed. 

The altitude curves for 0 to 40 kft were determined by plotting the altitude 
to sea level radiation ratios for data In Table 3.2 which were 0.1 Btu/ft 2 -sec 
or greater at any altitude. The result is shown In Fig. 5.1 with the solid 
curve selected to represent the current data, while the dashed curve shows the 1976 
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trend (Ref. 5). The curve of altitude ratios was used to determine points for 
construction of the family of altitude curves illustrated in Fig. 5.2. The 
higher curves generally apply to points with low heating rates which are shaded 
at oea level and get a better view of the plume as it expands. Points which have 
a good view of the plume at sea level are characterized by the lower curves. The 
lower curves decrease relatively slowly at low altitude to represent the per- 
sistence of the large Mach disk and shock structure in the plume. One of hese 
curves was assigned by the first altitude correction code letter to each point 
in the base region based upon (1) its view of the plume, and (2) the results 
reported for a similarly located point in Table 3.2. 

The second letter in the double letter altitude code specifies the radiation 
at 160 kft. It represents either a fraction of sea level heating or a fixed rate 
as shown at the top of Fig. 5.2. These are the same altitude correction codes as 
used in 1976 for the environment at 160 kft, and are determined using the vacuum 
plume predictions in Table 3.2. For points receiving significant radiation from 
the vacuum plume model, one of the ratios was used, while for points with 
negligible radiation one of the fixed rates was used depending on the predicted 
heating or the uncertainty due to recirculating gas and other undefined plume 
features. 

Although the SRBs are separated before 160 kft, an environment v/as specified 
at 160 k^t to allow interpolation between 40 kft and 160 kft. In assigning the 
160 kft correction factor for the SRB, the predicted heating of the inboard sur- 
faces of the SRB by the reversed gases from the SSME/SRB plume impingement was 
considered, and in some cases, this condition determines the environment specified 
at 160 kft rather than the anticipated radiation from the vacuum plume. Radiation 
from the SSME/SRB plume Intersection region was neglected. 
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[Altitude Corrections for 160 kft | 

Code Letter 
q i6ok /q SL 

q 16 ok( Btu/ft2 ' sec ^ 

JK L M N 0 P £ 

0.06 0.12 0.18 0.24 0.30 - 

0.1 0.2 
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39 .000 

3c9o.b 

1880. 5400 
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22.000 

19*2. 6 

2’. 010 

1976 . 0 

27.020 

3289.0 

2* .000 

3noi»? 

30.000 

3 yS 6 • 6 
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33 .000 

3196.1 
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3* .860 

7742.8 

4 7 .500 
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34 .9?J 

3524.5 

2116.2200 

.4601 

.4610 
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9*. *00 

399.2*0 

.4*50 

.7560 


94.719 

360.440 

.3600 

.7560 


94.729 

3056.860 

9.4350 

• 7560 


104. 105 

2958.360 

4.4160 

.7560 



41.670 96.010 


I%NE9 shock point TS 

• OOP 1608.700 


,000 1682 . 3*0 


1 * 9.3890 

170.5500 


20 


.7560 

.7560 
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FLOW FIELD iCCNT INUED I 


z 

4 1.670 


41.670 


41.670 


E T * 

98.010 


99.010 


180.000 



P 

T 

P 

FI 

10. 

roo 

1651.780 

117.0820 

.7560 

15. 

000 

1706.960 

126.4190 

.7560 

20. 

ooc 

1782.14Q 

15*. 3100 

.7560 

?5. 

OOP 

1890.590 

209.4470 

.7560 

30. 

000 

201 3.5?0 

283.9440 

.7560 

35. 

000 

2104.780 

354.1 550 

.7560 

40. 

000 

2050.840 

318.7250 

• 7360 

45. 

000 

1874.110 

215.3650 

.7560 

50. 

noo 

1627.490 

108.5350 

• 7560 

55. 

000 

1369.360 

51 ‘.4 810 

• 7560 

60. 

000 

1137.400 

25.9610 

.7*60 

65. 

non 

956.240 

13.2780 

.7560 

70. 

"00 

8o5.2»' 1 

9.4580 

• 7560 

75. 

000 

MG. 290 

3.9380 

.7560 

8". 

POO 

571.950 

1.5470 

.7560 

17?. 

944 

172.960 

.1560 

.7560 

17?. 

954 

2260.970 

.3010 

.7320 

173. 

8Q7 

2245.310 

.0010 

.7320 

Inner 

SHOCK POINT IS 

NUM4EP 18 


• 

000 

1608.700 

109.3890 

.7560 


000 

1681.050 

120.3510 

• 7560 

10. 

"00 

1682.000 

117.1570 

.7560 

i e . 

roo 

1706.330 

126.2040 

.7560 

20. 

oon 

1781.110 

1 56.8260 

.7560 

25. 

rco 

1889.360 

208,7010 

.7560 

30. 

POO 

2012.230 

283.19i0 

.7560 

35. 

noo 

2103.790 

353.3710 

• 7560 

40. 

roo 

2051.530 

319.1 710 

• 7560 

45. 

POO 

1336.200 

216.4260 

.7560 

5". 

000 

1630.350 

109.1880 

• 7560 

5*. 

roo 

1372.530 

51.8430 

• 7560 

60. 

non 

1141.060 

?6 .3 r 20 

• 7560 

t5. 

noo 

959.570 

13.4*40 

• 7560 

70. 

ooo 

867.0*0 

9.5020 

.7560 

75. 

noo 

711.730 

3.9580 

.7560 

80, 

ooo 

572.930 

1.5*30 

.7560 

• 

non 

1608.700 

109.3890 

• 7560 

c . 

ooo 

15° 1.310 

1 06.8610 

.7560 

1". 

ooo 

1694.890 

121.6220 

.7560 

15. 

POO 

167 3.6*0 

1 14.2730 

• 7560 

20. 

"00 

1733.220 

136.3160 

• 7560 

2 r • 

POO 

1831.310 

179.6830 

.7560 

30. 

"00 

1947.7*0 

244.3260 

• 7560 

35. 

POO 

2054.910 

314.9480 

.7560 

40. 

noo 

2097.760 

353.0090 

.7560 

45, 

"OP 

1973.270 

262.3*50 

.7560 

50. 

"CO 

1758.290 

151.2*30 

.7560 

55. 

oon 

1532.990 

43.4430 

.7560 

60. 

noo 

1317.610 

44.6540 

.7560 

65. 

noo 

1125.3*0 

23.3100 

.7560 
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LOW FIELD (CONTINUED) 


z 

ETA P 

T 

P 

FI 

4 1.670 

180*000 70. OOP 

968.990 

12.1230 

.7560 


75.000 

85M.820 

7.9390 

.7560 


60.000 

683.070 

3.2340 

• 7560 

_50.830 _ 

. . . .000 . .000 

17C3.0O0 

1 24 .4940 

.7560 


5.000 

1691.610 

120.5560 

• 7560 


lo.noo 

1697.690 

123.2650 

.7560 


15.000 

1766.510 

149.9090 

.7560 


20.000 

1872.280 

199.0690 

.7560 


25.000 

1986.930 

267.7930 

.7560 

. .. _ 

30.000 

2071.710 

327.6330 

.7560 


35. non 

2038. 5M0 

306.4720 

.7560 


MO. noc 

1883. M90 

209.0750 

.7560 


M 5, PqO 

1695.740 

123.7890 

.7560 


S'", noo 

1485.110 

69.2630 

• 7560 


55.000 

1273.550 

38.0540 

.7560 


60.000 

1081.620 

19.7610 

.7560 


62.557 

1056.240 

19.6880 

.7560 


62.567 

4446.670 

495.0000 

.7560 


7 C . 16? 

4446.670 

495.0000 

.7560 


inner shock point is 

NUMPEP 14 


50*830 

9.000 .000 

1703.000 

1 24.4940 

.7560 


r .Q00 

1691.710 

120.5890 

.7560 


10. non 

1697.290 

123.1180 

.7560 


l'.noo 

1765.870 

149.5870 

.7560 


20 . noo 

1871.560 

198.7410 

.7560 


25.000 

1986.040 

267.2760 

.7560 


3°. non 

2071.060 

327.1370 

.7560 


3*. noo 

2039.410 

307.1080 

.7560 


mo. non 

1884.610 

209.5210 

.7560 


M5. nuo 

169 7 s 3 T 3 

124.1990 

.7560 


50* 000 

1487.160 

69.7380 

.7560 


55. COO 

1275.900 

38.3210 

.7560 


57.M5M 

1176.510 

26.9290 

.7560 


57.M6M 

4446.670 

495.0000 

.7560 


67.72 3 

4446.670 

495.0000 

.7560 


inner shock point is 

NUMBER 13 


50*830 

9.500 .000 

1703.000 

1 24.4940 

.7560 


'.000 

169*. 720 

120.5930 

.7560 


i n *noo 

16°7 .250 

123.1010 

.7560 


l r> , noo 

1765.740 

149.5500 

.7560 


20.000 

1871 .440 

198.7030 

.7560 


25 . noo 

1985.990 

267.2170 

.7560 


30.000 

2070.940 

327.0800 

.7560 


35 . ''OO 

20 , 9.510 

307.1800 

.7560 


mo. noo 

1884.740 

209.5720 

.7560 


M 5.000 

1697.510 

124.2460 

.7560 


50.000 

1487.390 

69.7930 

• 7S6Q 


55.000 

1276.170 

38.3510 

.7560 


57.215 

11 85.440 

27.6090 

.7560 


57.225 

3293.900 

1 50.0930 

.7550 
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z 

50.830 

50.330 


50.330 


50.330 


50.330 


ETA 

R 

T 

P 

FI 

9.500 

69.492 

3291. 1*0 

1 33.2000 

.7550 


INNE 15 SHOCK POTNT IS 

NUMBER 13 


16.680 

.noo 

1703.000 

1 24.4940 

.7560 


5,000 

169 1.940 

120.6700 

.7560 


io. noo 

16«»6.340 

122.7630 

.7560 


15.000 

1764.150 

148.8090 

.7560 


20. OOP 

1869.830 

197.9480 

.7560 


25.000 

1934.050 

266.0300 

.7560 


3°. 000 

2069. 4R0 

325.9480 

.7560 


35.000 

2 041.520 

308.6410 

.7560 


40.000 

18 0 7 .4 n 0 

210.7240 

.7560 


4 r ' • 00 0 

1701.150 

1 25.1 850 

.7560 


50.000 

1492.090 

70.8820 

.7560 


54.939 

1305.270 

42.1080 

.7560 


54,449 

3282. 3H0 

172.34Q0 

.7550 


66.898 

3288.730 

1 64.5440 

• 7550 


Inner shock point is 

NUMBER 12 


24 • 3 10 

.000 

1703.000 

1 24.4940 

.7560 


5.000 

1692.320 

120.7970 

.7560 


10. ooo 

1694.830 

122.1483 

.7560 


15.000 

1761.510 

147.5730 

.7560 


20 . non 

1867.040 

146.6890 

.7560 


25. noo 

1930.840 

264.0610 

.7560 


30.000 

2067. 04J 

324.0880 

.7560 


3 r '. noo 

2044.850 

311.0590 

.7560 


40. noo 

18 9 2 .0 n 0 

212.8910 

.7560 


4 5. noo 

1707.140 

126.7370 

.7560 


50. noo 

1499.810 

72.6800 

.7560 


52.617 

13RC.370 

53.4020 

.7560 


52.627 

3265.920 

191.3680 

.7560 


65.414 

3268.740 

1 58.0440 

.7560 


inner shock °oint is 

NUMBER 1? 


31.960 

.000 

1703.000 

124.4440 

.7560 


5.000 

1692.820 

120.9700 

.7560 


10 . non 

1692.760 

121 .4260 

• 7560 


15.000 

1757.910 

1 45 .8900 

.7560 


20. noo 

1863.350 

194.9770 

.7560 


2 5. non 

1976.500 

261.3960 

.7560 


30. ngn 

2063.790 

321.6120 

.7560 


35. noo 

2049. 3TQ 

314.3270 

.7560 


40. noo 

1898.210 

215.8130 

.7560 


45.000 

1714.860 

129.0250 

.7560 


5''. non 

1510.120 

75.0970 

.7560 


51. *05 

14 3 5.170 

59.4220 

.7560 


51.815 

3250.410 

209.6840 

.7560 


6 r • 125 

325U.570 

2P8.6P80 

.7560 


inner shock rotnt is 

NUMBER 12 


45.550 

.000 

1703. OPO 

1 2 4 .4 R40 

.7560 


r .noo 

16«4.030 

121.3830 

.7560 


in. non 

1697.730 

119.5550 

.7560 



F L fy FIELD ( CONT INUEP ) 


z 

CT* 9 

T 

P 

FI 

50.833 

45.550 15.000 

1752.170 

143.6340 

.7560 


20.000 

1854.470 

1OQ.8720 

.7560 


25.000 

1966.210 

255.3710 

.7560 

, , . . „ „„ _ _ , 

_ 30.QQO 

2056.320 

315.9700 

.7560 


35.000 

2060.110 

372.0500 

.7560 


40.000 

1912.920 

222.7010 

• 7560 


45.000 

1730.670 

1 36.6800 

.7560 


50, OQO 

1531.4"0 

80.3890 

.7560 


54. *2? 

1341.140 

46.2380 

• 7560 


54. *37 

3254.840 

179.3030 

.7560 


6*. 565 

3255. 9?0 

1 71.2200 

.7560 


Innep shock POINT IS 

NUM4E° 1? 


50.830 

59.140 .000 

1703.000 

1 24.4 940 

.7560 


5.000 

1695 ,6 n 0 

121.9190 

• 7560 


lo.noo 

1681.160 

1 17,3750 

• 7560 


15.000 

1745.450 

141.1680 

.7560 


ZO. 000 

1842.820 

185.5190 

.7560 


25.000 

1953.070 

246.9590 

.7560 


30.000 

20*7.310 

308.9980 

.7560 


35. 000 

2073.950 

331.8900 

.7560 


4n.ooo 

1930.970 

231 .7480 

• 7560 


4". 000 

1750.600 

146.3570 

.7560 


50. non 

1555.440 

86.6360 

.7560 


55 . non 

1356.060 

48.2980 

.7560 


60.000 

1165.830 

25.5630 

.7560 


61.658 

1107.980 

20.9070 

.7560 


6 1.668 

3238 .830 

1 31.8760 

.7560 


74.526 

3225.4*0 

108.3590 

.7560 


IMNEP shock p otnt is 

NUMBEP 14 


50.93P 

72.7 30 .000 

1703. ono 

1 24 .4940 

.7560 


r .ooo 

1607.440 

122.5590 

.7560 


10.000 

1680.640 

1 16.8180 

.7560 


15.000 

1737.200 

1 78.1570 

.7560 


20.000 

1828.980 

1 79.1050 

.7560 


25. nQO 

1937.550 

277.4330 

.7560 


3°. non 

2037.570 

301 .4120 

• 7560 


35. TOO 

2090.260 

343.3460 

.7560 


40.000 

1947.660 

244.1160 

.7560 


45.000 

1774.050 

1 57.5130 

.7560 


5 r . non 

158 3.240 

93.8460 

.7560 


5 e • 000 

:38b. mo 

52.4770 

.7560 


6'"'. 000 

1194.750 

27.9880 

.7560 


6 r . non 

1022.750 

14.1710 

• 7560 


70. noo 

860.150 

7.3480 

.7560 


75.000 

701.160 

3.5110 

.7560 


7 r , 4 35 

637.380 

3.2740 

.7560 


7". 445 

31 3 7.620 

61 .0860 

• 7560 


84.61? 

3066.4OQ 

35.3300 

.7560 


Innep shock point is 

NUM8E9 17 


50.830 

"6.3 30 .noo 

1703.000 

1 24.4940 

.7560 
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(5 


FLOW FIELD (CONTINUED) 
7 ET ft 

50*830 86* 


50.830 99 • 9 20 


50*830 100*920 


I 5. 

non 

1699.610 

123.2870 

in. 

noo 

168 3.670 

117.8610 

15. 

OOP 

1727.960 

1 34.6880 

2D. 

noo 

181 3.480 

171 .8780 

2*. 

non 

1920.550 

226.9490 

3D. 

noo 

2027.740 

293.7350 

3*. 

noo 

2084.430 

338.6510 

40. 

noo 

1966.040 

257.5140 

4*. 

noo 

1799.470 

169.5210 

50 • 

oon 

1613.360 

101.6150 

55. 

noo 

1418.630 

57.0240 

60. 

non 

1226.080 

30.6630 

65. 

noo 

1056.520 

16.7250 

7H. 

000 

895.600 

8.8970 

75. 

000 

739.120 

4.4040 

- .._en. 

mo 

609.860 

2.1080 

85. 

POO 

514.130 

1 .1190 

9.1. 

ooo 

430.300 

.5850 

95. 

noo 

375.120 

.3530 

100. 

ooo 

329.350 

.2110 

104. 

690 

279.890 

.1200 

104. 

7Q0 

2785.890 

7.6020 

117. 

388 

2700.720 

5.8640 

inner 

SHOCK POINT IS 

NUMBER 

i • 

COO 

1703.000 

124.4940 

c . 

non 

1700.970 

123.7580 

10. 

noo 

1696.780 

1 18.9250 

15, 

non 

1717.810 

1 30.8940 

2°. 

POO 

1797.020 

1 64 .1640 

25. 

ooo 

1902. 89U 

216.0340 

30. 

ooo 

201 3.900 

284.3610 

36. 

noo 

2071. 1*»0 

328.5240 

40. 

noo 

1985.160 

271.2110 

46. 

noo 

1825.590 

181.7100 

5 

non 

1644.160 

109.5110 

56. 

noo 

145 1.860 

61 .6920 

60. 

POO 

1260. 5T0 

34.7830 

66. 

ooo 

10° 1 .230 

19.5690 

70. 

POO 

931 .660 

10.5370 

75. 

oon 

777.710 

5.3700 

80. 

POO 

64 1 .820 

2.6150 

85. 

oon 

544.580 

1.4330 

90. 

noo 

459.880 

.8160 

188. 

217 

141.500 

.0120 

188. 

227 

2262.510 

.0010 

1 89. 

535 

2256.880 

• 0010 

Inner 

SHOCK POINT IS 

NUMBEP ; 

• 

OOO 

1703.000 

124.4940 

% 

POO 

1701.020 

123.7730 

10. 

noo 

1687.010 

1 19.0030 


27 


FI 

• 7560 
.7560 
.7560 
.7560 
.7560 
.7560 

• 7560 
.7560 
.7560 
.7560 

• 7560 

• 7560 

• 7560 
.7560 
.7560 
.7560 

• 7560 
.7560 
.7560 

• 7560 
.7560 
.7550 
.7530 

.7560 

• 7560 
.7560 
.7560 
.7560 

• 7560 
.7560 

• 7560 
.7560 

• 7560 
.7560 

• 7560 
.7560 
.7560 
.7560 
.7560 
.7560 
.7560 
.7560 

• 7560 
.7370 
.7370 

.7560 

.7560 

.7560 
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CW FIELD 
2 

50.830 


SO* 830 


60.030 


6 0.000 


ETA 

100*920 


180.000 


,0 00 


11.000 


:oi 




.... 


p 

T 

P 

FI 

is. 

000 

1717.050 

1 30.6070 

• 7560 

?o. 

000 

1795.790 

1 63.5880 

.7560 

2*5. 

000 

1901.590 

215.2300 

• 7560 

-30. 

000. 

.2012.530 

283.5180 

.7560 _ 

38. 

000 

2070.260 

327.8170 

• 7560 

40. 

000 

1986.560 

272.2100 

.7560 

4S. 

000 

1827.500 

1 82.5950 

• 7S60 

5 r . 

000 

1696.910 

1 10.0850 

• 7560 

55. 

000 

1959.280 

62.0330 

.7560 

60. 

000 

1263.200 

35.1470 

.7560 

6S. 

000 

1093. 7f0 

19.7790 

.7560 

70. 

000 

939.270 

10.6590 

.7560 

75. 

000 

780.510 

5.4420 

• 7560 

80 . 

000 

699.n0 

2.6540 

• 7560 

85. 

000 

596.810 

1.4570 

.7560 

90. 

OQQ 

__ 96 2.180 

• 8360 

• 75.60 .. 

. 

Ooo 

1703.000 

124.4940 

• 7S60 

5. 

noo 

1709.950 

125.0280 

• 7560 

10. 

000 

1698.500 

122.9040 

• 7560 

15. 

ooo 

1689.900 

118.1760 

.7560 

20. 

Ooo 

1793.890 

140.9490 

• 7560 

25. 

ooo 

1838,920 

183.8410 

♦7560. 

3 n . 

ooo 

1995.620 

242.1980 

.7560 

35. 

coo 

2033.920 

299.0950 

.7560 

90. 

ooo 

2056.310 

320.1230 

.7560 

9 5. 

000 

191 5.670 

226.3790 

.7560 

5". 

ooo 

1748.930 

1 46.9170 

.7560 

5 5. 

ooo 

1570.000 

89.7860 

.7560 

60. 

non 

1388.200 

52.6610 

.7560 

65. 

coo 

1214.080 

30.1430 

.7560 

70. 

ooo 

1058.320 

16.8520 

.7560 

75. 

ooo 

91 3.660 

9.7680 

.7560 

80 . 

noo 

774.100 

5.5490 

• 7560 

88. 

000 

663.670 

3,2100 

.7560 

9 n , 

non 

574.400 

1.9460 

.7560 

. 

ooo 

1689.680 

119.8470 

.7560 

5. 

non 

1698.283 

123.1150 

.7560 

10. 

ooo 

1737.270 

1 38.8200 

.7560 

15. 

ooo 

1816.100 

1 72.9550 

.7560 

20. 

000 

1915.400 

222.7090 

.7560 

28. 

noo 

2008.320 

280.9250 

.7560 

30. 

ooo 

2028.340 

298.4010 

.7560 

38. 

roo 

1950.710 

251.3390 

.7560 

9°. 

ooo 

1806.850 

173.7940 

.7560 

95. 

ooo 

1646.150 

108.9090 

• 7560 

50, 

,000 

1470.300 

64.3510 

.7560 

55. 

non 

1292.070 

38.9160 

• 7560 

5<5. 

982 

1129.780 

22.7840 

.7560 

5^>. 

992 

4320.000 

540.0000 

.7560 

7* . 

162 

432C.000 

540.0000 

.7560 

nep 

SHOCK POINT IS 

NUMBEP 13 


. 

noo 

16B9.68Q 

119.8470 

.7560 
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rtcw FIELD (CONTINUED) 
l ETA 

60.000 2b« 290 


6 0. GOO 32.440 


i ' 


i - 


J 

«► 

I 

I 

I 


60.000 46.220 


6 0. 0 JO 6C.0 2Q 


30. 

nor 

2026 

.990 

297 

.0750 

3 6. 

OOD 

1958 

.090 

256 

.0820 

40. 

non 

1818 

.540 

1 78 

.0190 

.4 5. 

000 

1657 

.590 

111 

.7530 

4*». 

927 

1485 

.340 

67 

.9760 

40. 

937 

3018 

.750 

22 7 

.8930 

65. 

21? 

3021 

.260 

224 

.6840 

inneo 

5H0CK POINT IS 

NUMPEP 

• 

000 

16 8 9 

.680 
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